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\ We study the ground-state properties of a quasi-two-dimensional Bose- 

\ Einstein condensate of indirect excitons, which are confined in an anisotropic 

\ harmonic potential. Incorporating the interactions, we calculate the order 

' parameter variationally. The difficulties in the detection of a Bose-Einstein 

I condensate are also discussed, along with possible ways which would overcome 
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them. 



_ I. INTRODUCTION 

B. 

^ , Excitons are bound states of electrons and holes which form in excited semiconductors 

c : 

O . [1], very much like positronium atoms. Since excitons consist of two fermions, at least as 

O 

long as the mean spacing between them is much larger than their Bohr radius, they are 
^ ! expected to behave as bosons [2]. While the phase transition to a Bose-Einstein condensed 
phase has been realized experimentally in vapours of alkali-metal atoms [3], the formation 
of an exciton Bose-Einstein condensate (BEG) has turned out to be more tricky. A lot of 
experimental effort has been made on the formation of a BEG of excitons in GU2O [4,5], 
and in quantum wells [6]. Excitons resemble atoms, with one major difference being their 
mass, since they are much lighter particles and as a result they exhibit quantum-mechanical 
effects more easily. In that respect, excitons are more advantageous than atoms for forming 
a BEG. On the other hand, their basic disadvantage is their finite lifetime due to various 
decay mechanisms [7,8], which destroy them. 

For a gas of bosons of mass m and density n, the phase transition to a Bose-Einstein 



condensed phase occurs when their thermal de-Broghe wavelength becomes comparable to 
their interparticle distance, which imphes that ksT^ ~ U^ri^/^ /m, where is the critical 
temperature. For Tc being on the order of liquid-helium temperature, i.e., 1 K, and for 
m being on the order of the electron mass, n turns out to be on the order of 10^^ cm~^. 
For such densities, the boson-boson spacing is ~ 500 A, while typically the exciton Bohr 
radius is smaller than this length scale, and therefore the excitons are expected to behave 
as point-like (bosonic) particles. 

Recently in an interesting paper Butov et al. [9] (see also [10]) reported the formation 
of a degenerate gas of excitons in a quasi-two-dimensional geometry in coupled quantum 
wells. More precisely, in the experiment of Ref. [9] indirect excitons formed between two 
parallel layers, after an electric field was applied perpendicular to the two planes. Then, 
spatially resolving the photoluminescence that comes from the layers, isolated bright local- 
ized spots were observed. These spots were attributed to the recombining excitons which 
have concentrated in regions of local minima (lakes) that exist due to irregularities in the 
heterostructures. The presence of these local minima is important, since the excitons con- 
centrate in them. Furthermore, the exciton lifetime is long due to the small overlap between 
the electron and the hole wavef unctions, and also the coohng due to the scattering of excitons 
with acoustic phonons is enhanced because of the relaxation of the momentum conservation 
along the direction of confinement. All the above effects favor the formation of a cold and 
dense gas of excitons, and indeed signs of degenerate behaviour in these lakes were observed. 

Motivated by this experiment, we examine in the present paper the ground-state prop- 
erties of a Bose- Einstein condensate of excitons in such a system [11]. Our goal is to demon- 
strate some important physics which, in addition to other things, would allow the detection 
of a Bose-Einstein condensate and not to give a detailed quantitative description. On the 
other hand, our results should be reliable at least as order-of-magnitude estimates. In that 
spirit, and since the real system is quite complicated, we make some simplifying assump- 
tions. For example, we neglect the finite lifetime of excitons, since typically this is larger 
by one or two orders of magnitude [12] than the thermalization timescale [13]. In addition, 
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the rate of non-radiative, Auger-like processes is expected to be suppressed because of the 
spatial separation between the electron and the hole of each individual exciton [14]. 

In what follows we examine in Sec. II the interactions between the excitons and in Sec. Ill 
how the statistics due to the Pauh principle between the electrons and/or the holes of 
different excitons affect them. We then make some estimates in Sec. IV, where we model the 
in-plane confinement with a harmonic potential. Using a variational ansatz for the order 
parameter we calculate in Sec.V the energy and the width of the cloud in the case of an 
asymmetric confining potential. Finally we discuss in Sec. VI how one could get evidence 
for the existence of a BEC, while Sec. VII summarizes our results. As shown, the width of 
the condensate is expected to be of the same order as that of the thermal cloud and for 
finite temperatures, that would obscure the detection of the condensate. In that respect, 
an important point of our study is that (similarly to the atomic condensates [15]) as long 
as the trapping potential is anisotropic, the kinetic-energy distribution of an expanding 
BEC (upon release of the excitons from the trap) is also anisotropic, while the kinetic- 
energy distribution of a thermal gas is always isotropic. Therefore, any experiment which is 
sensitive to the anisotropy of the kinetic-energy distribution of the excitons could be used 
as a "smoking gun" for the detection of a BEC. 

II. INTERACTION BETWEEN THE EXCITONS 

As mentioned above, the interactions play a very important role in this system. Numer- 
ous studies have examined this problem - mostly the effect of the reduced dimensionality 
and of the band structure of specific materials. It is not the purpose of the present paper to 
investigate such questions in detail, but rather to follow a method that captures the essential 
physics we want to demonstrate. 

It is important to remember that the electrons and the holes in this system are spatially 
separated and as a result the corresponding electron-hole pairs (that form in the limit of 
dilute densities that we consider) have a dipole moment. As shown in Ref. [16] in such 
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a system the effective exciton-exciton interaction is dominated by tfie one tfiat reduces to 
tfie usual dipole-dipole term for large exciton separations and this is the approach that we 
follow here. Viewed in another way, this approach is phenomenological. While it is certainly 
reliable for large separations between the excitons, it is less rehable for short distances, 
where one needs to worry about the exchange terms between the electrons and the holes of 
different excitons. However, the dipole moment of the excitons dominates the interaction 
and our approach is certainly valid for low exciton densities; a similar approach was followed 
in Ref. [17]. References [18] have studied a problem close to the present one in the context 
of atomic condensates which are subjected to an external electric field. 

More precisely, let us consider the two-exciton problem, i.e., two electrons confined on 
a plane and two holes confined on a parallel plane, with the planes being separated by a 
distance D. As in the experiment of Ref. [9] we assume that D is on the order of the Bohr 
radius of a single exciton (this is a reasonable thing to do, since if D ~ the electron- 
hole bound state is essentially a p-like hydrogenic state). Fixing the position of, say, the 
two electrons at an infinite distance, in the lowest state each hole binds with one electron 
and two excitons form, with the total energy of the system being twice the binding energy 
of a single exciton. If the two excitons are brought to a finite distance from each other, the 
interaction between them is given by 

e2 / 2 1 1 \ 

^^^-^^^ = 7(,F^- |r-r- + D| - [r-r--D[ J' 

where e is the dielectric constant, and r — r' is the vector connecting the electrons or the holes 
of the two excitons. In our study we assume for simplicity that the exciton dipole moment 
is parallel to the electric field (i.e., perpendicular to the x-y plane) and we neglect the 
fluctuations around this orientation. For small distances between the two excitons, |r — r'| -C 
Eq. (1) reduces to the usual Coulom interaction between the two electrons and the two 
holes, V{v — r') = 2e^/ (e|r — r'|). For large distances between the two excitons, |r — r'| » D, 
Eq. (1) reduces to the usual dipole-dipole interaction, V{t—y') = e^Z^^(l— 3 cos^ ^)/^k ~ '^1'^' 
where 9 is the angle between D and r — r'. In our case, since the system is quasi-two- 
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dimensional, Q 'Kjl. Therefore V{x — r') is purely repulsive, which guarantees that there 
are no instabilities against the formation of e.g., molecules. 

III. EFFECT OF THE PAULI PRINCIPLE 

When the separation between the excitons becomes of order a^, one has to worry about 
the Pauli principle between the electrons and/or the holes. If the spins of either the elec- 
trons, or the holes are parallel, the relative wavefunction between the two excitons is highly 
suppressed for values of |r — r'| smaller than as (on the other hand, this effect is absent 
when both the electrons and the holes have antiparallel spins). As a result, the excitons are 
not allowed to get to distances shorter than |r — r'| ~ as- Therefore, for exciton-exciton 
separations on the order of a^, the true many-body wavefunction behaves in a similar way 
(since for the low densities we consider, the effect of other excitons is negligible) [19]; how- 
ever our mean-field approach ignores these correlations at such small scales. For this reason, 
to calculate the interaction energy when the Pauli principle is active, we will use Eq. (1) 
assuming that the two excitons cannot get to distances closer than a^. It is important to 
point out two things concerning these implications that result from the Pauli principle: (i) 
as shown in detail below, the energy of the system is not very sensitive to the choice of the 
cutoff length (which we choose to be for simplicity), and (ii) the bosonic nature of the 
electron-hole pairs is not affected, as long as the mean exciton-exciton spacing is larger than 
Ob, which is indeed the case in the experiment of Ref. [9]. 

IV. ESTIMATES 

Before we calculate the order parameter of the condensate, it is instructive to make some 
estimates. Let us assume that the excitons are confined on the x-y plane within a lake of 
radius d. For low enough densities the dipole expansion can be used and in this case the 
energy per particle of a BEC of excitons is (apart from numerical factors of order unity). 
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where N is the number of excitons. In the above equation we have assumed that the excitons 
are trapped in a harmonic potential of frequency uj (corresponding to the confinement in 
the lake). Equation (2) can also be written as 

where = j y.^^ is the exciton Bohr radius, with [i^ — memh/m being the reduced 
mass of the effective electron mass me and the effective hole mass m/j, and m — rrie + nih- 
For ~ O.lmo, where ttiq is the electron mass and e ~ 10, then as ~ 100 A, which 
is ~ L>. Under typical conditions, N ~ 10^, the first term on the right side of Eq. (3) 
is negligible. Under this assumption, the confinement due to the lake is balanced by the 
(repulsive) interactions between the dipoles, while the kinetic energy (due to the Heisenberg 
uncertainty principle) is negligible, very much like the Thomas-Fermi approximation in 
atomic condensates [20], 

and thus the radius of the condensate is 

do/ao,c^ci{ND/aBy/\ (5) 

with Cose {h/mujY^'^ being the oscillator length and ci being a constant of order unity. 
Also, 

S{do)/nuj = C2{ND/aBY'\ (6) 

where C2 is a constant of order unity. A useful formula is the one that compares the radius 
of the condensate do with the radius of the thermal cloud R. Since muj'^B? ~ ksT, and 
assuming that ksTc ~ hcuN^^'^ [3], 



Remarkably, the ratio d^/ R does not depend on the number of excitons N . Turning to the 
speed of sound c, mc^ ~ (y), where {V) ~ Ne^D/ed^ is the interaction energy per particle, 
and thus c ~ {Uj mdQ){N D j asY^^ ■ Also, the coherence or healing length ^, satisfies the 
equation /2m^ ~ (V), and thus ^/rfo ~ {aB/NDy^'^. Finally the critical frequency Qc for 
creating a vortex state is ~ (h/mdl) \n{do/^). 

Let us make some numerical estimates now. Since the radius i? of a thermal cloud of 
excitons in the lakes is 10 //m for T 1 K [9], the formula mu'^E?' ~ ksT imphes that 
uj ~ 10^ Hz for m being on the order of O.lmo- Given this value for to, the oscillator length 
Cose is ~ 1 /xm, and /io; ~ 5 x 10""^ meV. Therefore, since D ~ a^, then ND/gb ~ iV ~ 10^, 
and the size of the condensate do is ~ 30 //m, the surface density n2D is ~ 3 x 10^° cm~^ 
and the mean exciton-exciton spacing ~ d^/N^/'^ is ~ 300 A. The fact that n2D (i% is ~ 0.03, 
i.e., much less than unity, where n2D is the surface density, justifies the treatment of the 
excitons as bosons and the use of Eq. (1) as the appropriate expression for the interaction 
between the excitons. Turning to the energy per particle of the condensate, this is ~ 1 meV. 
The transition temperature Tc is ~ 5 K, while for T = Tc, Eq. (7) implies that do ^ R. Also 
the speed of sound c is ~ 10^ cm/s, and the coherence length ^ is ~ 100 A, which gives 
the typical size of vortices [10]. Since ^/c?o ~ 10~^, the superfiuid properties of this system 
should resemble those of helium (and not those of nuclei). Finally the critical frequency Qc 
for creating a vortex state is ~ 5 x 10^ Hz. 

V. VARIATIONAL CALCULATION OF THE ORDER PARAMETER 

Let us now get more precise. If ^'(r) is the order parameter of a BEG of excitons, as 
long as the estimates made earlier are accurate, it is enough to consider only the relevant 
terms in the Hamiltonian, 



where Vt(r) is the trapping potential and is the chemical potential. We will now use the 
variational approach to calculate the order parameter and the energy. For simplicity we 
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restrict ourselves to two dimensions [21], i.e., on the x-y plane and assume that the trapping 
potential has the form 

Vt{v)^]^muj\x' + W). (9) 
with uj — ujxi and A = LOyjuj^-, and also 

= ^1/2^. r, u/2 ^^P(-^V2a^ - 1/V2a5), (10) 

where and Uy are variational parameters. Then, the energy per particle is, 

= ^ (/ Vt{r)\^{r)\^d\ r')|^(r)n^(r')|'c^VdV) . (11) 

Starting with the expectation value of T4, /N = ma;^(a^ + A^a^) /4. Turning to the 

interaction energy 

6 /'it/\i9it//\i0( 2 1 1 



(^) = ^ / l^(r)n^(r')r - 1 — - 1 = r d\d\\ (12) 

^ ' 2NeJ ^ ^ ^' ' ^ ^' V|r-r'| |r-r' + D| |r-r'-D|y ' ^ ^ 

where, as we argued earlier, |r — r'| > gb- To evaluate the above integral, we use the relative 
coordinates ri2 = r — r', and R12 = (r + r')/2. Equation (12) then involves an integration 
over R12 and ri2, 

(V) = , f / e-^«?^/(^) d'Rr, f-(l- , \ e-^-^(^')/^ dV,2, (13) 

where /(0) = (cos0/a^)^ + (sin^/a^^)^. Since D -C ax,ay, we take the exponential in the 
integral over ri2 to be equal to unity, and thus 

Jas ri2 \^ ^rl2 + Dy ri2 \^ ^rj^ + Dy 



^2TT{^Jal + D^-aB). (14) 
Combining the equations above, the energy per particle is 



expressed in terms of the two variational parameters ax, ay, and A. For A = 1, then ax — 
ay = d, f = the integral in Eq. (15) is equal to t^cP /2, and Eq. (15) gets the simple 

form 



1 , Ne\^al + D--as) 

8 = -mw d H — . (16) 

z ed'^ 



As Eq. (16) implies, the energy of the system scales as {\ja\ + D'^ — a^)^/^, which shows 
that indeed £ is not very sensitive to the cutoff length for ~ D. Furthermore, if the 
cutoff length in Eq. (14) is set equal to zero, 

£ = -mu d + (17) 

in agreement with the estimates of Eq. (4). Equation (17) implies that Ci in Eq. (5) is equal 
to {2m/ jixY^^ , and C2 in Eq. (6) is equal to {2m/ jixY^"^ , where m — me-\- m^- 

Returning to the more general problem of A ^ 1, we have minimized the energy of 
Eq. (15), setting for simplicity the cutoff length equal to zero. Measuring ax and ay in the 
second term of this expression in units of Oosc, the energy scale 2Ne'^D /i^ea^^ can be written 
as Tujj{2ND /T^aB){m/ jjLx)- We choose {2N D / T^aB){m/ jix) to be equal to 0.5 x 10^. We also 
consider the value A = 0.8, which corresponds to a trapping potential with equipotential 
lines which are elongated along the y axis, with a ratio of the oscillator lengths aosc,x / (iosc,y 
equal to A^/^ ^ 0.894. Minimizing numerically the energy, we find that £^/^a; 1.119 x 10^ 
with ax/aosc ~ 33.3 and ay/a^sc ~ 42.0. It is interesting that the ratio a^/ay ^ 0.793 is 
smaller than A^/^. Therefore, the interactions enhance the anisotropy of the cloud due to the 
trapping potential, (i.e., they make the interacting cloud more elongated along the y axis), 
as compared to a non-interacting gas, which would result from the uncertainty principle. A 
similar effect has been observed in atomic condensates with contact interactions [15,20]. 

VI. DISCUSSION OF THE RESULTS AND EXPERIMENTAL RELEVANCE 



In the last part of this study we examine how the above results can be used in practice. 
Getting indisputable evidence and clear signatures for the existence of a BEG of excitons is 
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not at all a trivial problem [23] , since even for small but finite temperatures the condensate 
co-exists with the thermal cloud. More precisely, as Eq. (7) implies, the width of the con- 
densate is, even for T = Tc, comparable or even larger than that of the thermal component, 
and this ratio gets even larger as the temperature decreases. This is opposite to the atomic 
condensates, where typically the thermal cloud is much broader than the condensate. This 
fact would certainly obscure the detection of a BEC of excitons in the present system. 

In the alkali-metal atoms, in the first experiment where the formation of a BEC was 
reported [15], the atoms were confined in a disk-shaped trap. Following their release from 
the trap, for very low - essentially zero - temperatures, the kinetic-energy distribution of the 
atoms was observed to be anisotropic, while for higher temperatures (but still lower than 
Tc) a more broad isotropic background was also observed, corresponding to the thermal 
component. Above T^, only the isotropic component remained. Therefore, the anisotropic 
expansion of the cloud is a unique characteristic of a BEC, and does not show up in a thermal 
cloud. 

Turning to the present problem, ideally one could think about releasing somehow the 
excitons from the lake and watching them as they expand. Actually, such a possibility 
has already been studied theoretically [24]. In such an experiment, any sign of anisotropic 
expansion of the cloud on the x-y plane below some temperature (the transition temperature 
Tc) would provide strong evidence in favour of the presence of a BEC of excitons. On 
the other hand, from the estimates made earlier, one can see that the typical velocity of 
expansion is expected to be ~ (£^/m)^/^ ~ 10^ cm/s. Since excitons with velocities up to 
^ 1.4 X 10^ are radiatively active [25], the expanding excitons may or may not be radiatively 
active. Additionally, since the velocity of expansion is close to the sound velocity of the 
material (~ 3.7 x 10^ cm/s), the interaction with phonons will slow them down. Apart from 
these possible difficulties, more generally one should remember that any experiment that is 
sensitive to the directionality of the kinetic-energy distribution of the excitons could be used 
to provide evidence for the formation of a BEC. 
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VII. SUMMARY 



To summarize, in this study we examined the ground-state properties of a quasi- 
two-dimensional Bose-Einstein condensate of indirect excitons, which are confined in an 
anisotropic harmonic potential. Incorporating the interactions we made estimates for the 
size of the cloud, the energy per particle, the density, the speed of sound, the coherence 
length, and the critical frequency for creating a vortex state. We also calculated the order 
parameter variationally within the mean-field approximation and we examined the effect of 
the Pauh principle. Concerning the detection of a BEC, since the size of the condensate is 
expected to be of the same order as that of the thermal cloud, we examined the possibil- 
ity of releasing the excitons from an anisotropic trap, which would result in an anisotropic 
kinetic-energy distribution and it would provide evidence for the presence of a BEC. An 
interesting result is that the interacting gas exhibits this anisotropy in a more pronounced 
way as compared to an ideal gas. 
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